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A77-S1267 i‘ Numerical Laplace back-transformation for 

determining system response: of linear systems in the time domain 
(Numerische Laplace-Riicktransformation zur Bestimmung dcr Sys- 
temantwort linearer Systeme im Zeitbcreich). R. Friedrich and VV. 
Orevvelovv (Rostock, Universitat, Rostock, East Germany). Rostock, 
WUhehn-Pieck-Universitat, Wissenschaftliche Zeitschrift, Mathe-' 
maiisch-Naturwisscnschaftliche Reihc, vol. 25. no. 8, 1976, p. 
887-889. 7 refs. In German. 

This paper present: a program (or the numerical Laplace return 
transformation. The algorithm described is based on the method of 
breaking the picture function down to form partial fractions which 
are then transformed back separately. The sum of the resulting 
partial functions is the sought time function. Any problems caused 
by the forms of a system of equations being unfavorable due to the 
computer configuration are restricted to a large extent by appropri- 
ate standardization using an auxiliary parameter. The practical limits 
for use of the program are reached when the degree of the 
denominator of the picture function lies within the range of seven to 
eight. (Author) 



Thu Abctraiet l^uBllshwl In 
ilOfEfUIAriONAtrAEfK^ACELABSTRACTS - 


Issu^ by 

TECHNICAL INFORMATION SERVICE 
AMERICAN INSTITUTE OF AERONAUTICS AND ASTRONAUTICS, INC. 
750 THIRD AVENUE, NEW YORK, N.Y. 10017 


SYS 

MAI 

APR I 




Numerischc Laplacc-Rucktransformation zur Bcstimmung dcr 
Syslemanlwort lincarer Syslcme im Zeitbereich 


von Half Frioclrich und WolfRanH Oi«-vvflt>w 


1. Einleituni; 

KinsdiwintiunKsvorKanKe spii-len in tier Tcchnik cine 
wichtiKe Rollc. Die mathematische Heschrcibunc clieser 
VorRiince crfolgt durch DiffcrentialRlcichunKon. Eines 
dcr wirksamstcn Hilfsmiltel fur den Technikcr bci dcr 
LbsunR dieser AufRaben ist die Laplace-Transformation. 
Die zuRchiiriRe Laplacc-Rucklransformation bereilct 
jcdudi oft wcRen dcs hohcn Rcchcnaufwandes Schwie- 
riRkcitcn. Insbesondere <Iann, wcnn dicsc EinschwinR- 
vorRiinRc fiir cine RrofJcre Anzahl vcrschiedener 
Zcitwcrtc borechnet wcrdcn sfillen. crReben sich in 
pruktisch wichtiRen Fiillcn umfanRiciche Rechenarbei- 
tcn. 

Aus diesem Grunde schicn cs sinnvoU, cin ProRiamm 
fur cincn DiRitaIrcchncr zu cntwickeln. mil dcm es 
miiRlich ist, au.s eincr Bildfunktion die zuRchiiriRc 
Funktion im Zeitbereich zu bestimmen und cine Werte- 
talx'lle, cnlsprechend vorRcRcbener ForrlerunRcn, aus- 
di'ucken zu lassen. fm VordeiRi uncl Manden die wich- 
tiRen aperiodischen .Syslemreaklioncn ObetRanRsfunk- 
lion, Gcwichlsfunklion und Rampcnantwoiifunklion. 
Die VerwendunR harmonist+ier SiRnale .sollle <d)enf.ills 
miiRlich scin. 

Ein wichliRes AnW' ndunRsRcbie'. fiir tiieses ProRiamm 
ist die Syslemi<ientifikali<>n, rl;; mil (lessen Ililfe die 
Giite verschiedener Verfaiiren im Zeitbereich Rut iiber- 
pruft werden kann, Des ueiUren kann das ProRramm 
beim Unlersuchen des KuhruiiRs- und Sliirverhallens 
Rcsschlossener linrsaier HeRclkreise ebensn erfolRieich 
einRcsclzl werden wie beim Beslirnmen der Ein- 
schwinRvorRiinRc in elekironischen SthallunRen. 

Zur IzisunR stand der Kicinrechner Cellalrnn C 1120S 
zur VerfufiunR. 


2. OrundlaRen 

Ein linenres. System mil konzentrierten Par.'imetern 
und dem EinRaiiRssiRnal x,, t bzw. dcm AusRanRssiRnal 
X., I wild durch die DifferenlialRleichunR 

dx,;’ dx.," ' 


dx m 


+ b.„ 


dxj' 
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beschrieben. Fiir vcrschwindende AnfanRsbedinRunRcn 
liefert der Differenliatinnssalz der Laplace-Transforma- 
tinn 


^ b,nP"' + »’n 
P :i„ p" f a„ 


,l>"’ '4 ... i b„ 
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G(p) ist dabei die UbertraRunRsfunklion des linearen 
Systems in Polynomform. Fur die BelanRc der Rcrc- 
lunRstechnik ist jedoch die soRcnannte Zeitkonslanlen- 
form der ObertraRupRsfunktion besser ReeiRnet, Be- 
schriinkt man sich aus Ober.sic+illichkeilsRriinden vnr- 
ersl auf rccilc Pole und .Vullslellen der rtitionalen 
Funkticjn nach Cl. (2), so folRt 


G'p =K, o'- ' 
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Im Fall konjuRierl knmplexer Pole Oder Nullslcllen 
werden jeweils zwei Klammeraii.sdrucke 
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zusammcnRelabt. Unter Hmzunahme der I.aplaceliiins- 
formation des EinRnnRssiRntils crhiilt man 

X,, p, = X^. p G |) . 

Zur Lapl.ace-Rucktransrorm.it ion von X.^ p kann man 
den Residuensatz heranziehen, jedoch miissen dann fiir 
mehrere Pole Differenlialionen durchRefiihrl werden. 
die auf dem DiRilalrechner .sehlechl zu realisieren sind. 
Deshalb ist die ZcrlcRunR der Funklion X., (p) mich 
G1 (■'{)— (~i) in eine .Summe von Parlialbriichen das bes- 
scr ReeiRnete Verfahren. 

Man erhall so Briiche der Form 
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Dii-se nriii'hc sind tcicht ruc-ktransformicrbar f2l. die 
Sumnie diT sn enlstaudenen /oiU'iink)i<inen ist die Ke- 
sudilt- SysIrmantuYdl im Zcillwreif-h. 

AIs AufiialH' vcibleibl. die Konstanten, cleren Anzahl 
Uleij+i dem NenneiKi'ad von p — im folRcndcn 
urad — isl. zu lieslimmcn. 

Dazu stehen die Melhfide des KoeffizicntenvcrKleichs. 
tlie CirenzwcrtmethtKie und die EinsetzunRsmethode zur 
V’erfuRunR. Die EinsetzunRsmethode ist dabei am leich- 
li-stcn proRrammierbar und wurde deshalb den anderen 
Verfahren vorRezoRen. Man <Thiilt (lurch das Einselzen 
von Read verschiedene p-U'ei ten Rrad linear unabhiin- 
RiRe CleichunRen zur IlerwhnunR der Read Konstanten. 
Die ParlialbruchzerlrRunR in VerbindunR mit dem Ein- 
setzunRsveifidiren als einer ret+meri.st!ien I.osunRS- 
rnethiide bildet das KernstCiek des im I’roRiamm ver- 
wendeten IxisunRswcRes. 


ProRrammbeschreibung 

Nach der EinRabe der erforderlichen aktuelien Parameter 
(Ol)ertraRunRsfunktion. Resuchte Systemreaktion. An- 
forderunRcn an die auszudruckende Werlelabelle) wird 
im Hechner auf Giund dierser EinRabewerIc die allRe- 
meine AnsatzRleichunc zur ParlialbruchentwicklunR fie- 
l)ildet. In dieser CleichunR. auf deren rechten Seitc die 
zu ermittelnden Koeffizienten der Partialbruche stehen. 
werdon nadi einer bestimmten BildunRsvorschrift Rrad 
verschiedeiH! Werte fiir p einResetzt. Auf die.se Weise 
erhiilt man Rrad linear unabhnnRiRe CleichunRen zur 
IJestimmunR der Rrad unbekannten Konstanten. Diese 
(j'leiehunRen bilden den AusRanRspunkt eines Unter- 
proRi'amms zur laisunR eines linearen GleichunRs- 
systems. am Er.de dieses UnterproRramms stehen die 
Resuchten Konstanten. Diese Konstanten werden nidit 
erst wieder d('n entsprechenden Partialbruchen. sondern 
Rleieh den I.aplaceruektransformierten dieser Partial- 
tniiche zuReordnel. Dadurch erhiilt man Rrad einzelne 
Z(ritfunktionen; die .Summe der Zeitfunk'ionen ist dio 
Resuchte l.aplaceruektransformier..! von X„ (p). die 
.SystemauSRanRsfunklion im Xeitbereich. Die.se Funkti.m 
X., (t) wil'd au.sRedruekt. ansehliclSend werden entspre- 
diend den eiiiReRebenen VVerten fiir die Schrillweite 
und obere HeiechnunRSRi enze Zeilwerlc und die dazu- 
RehbriRen Fiinktionsvverte erniitlell und au.sRcdruekl. 

Dieser prmzipiell riehtiRe I^isunRswcR weist jedoch 
bei der praktisdien RealisieruriR im Kochner we.sent- 
liche MiinRcl beziiRlich der GenauiRkeil der zu bestim- 
menden Konsl.mten und damit der GenauiRkeit der zu 
bestiniirenden /.eitfunklionim auf. Die.se Fehler sind 
;iuf die PildunR i eehenterhnisch unRiinstiRer Koeffizien- 
ienmalrizen. die den AusRanRspunkt des Unlerpro- 
Riamms zur Iz.sunR des linearen GleichunRs.systems 
bilden. zuruekzuftihren. Bei BenuIzuaR der Ein.selzunRs- 
melluide mid bei einer im Prinzip willkurli(k fest- 
Relt'Rlen BildunRsvorschrift fiir p erhidt man hiiufiR an- 
niiliernd linear .djhiinc.iRe GleiehunRe.n. In dicsem Fall 
befindet sieh der Heehher im UnlerproRr.'imm zur Izi- 
siiUR eines linearen GlcichunR.ssystems in einem relativ 
labileii Bereich. d. h. es kiinnen sehr R.-ol)e Reehenfeh- 
ler auftreten. Ilinzu kommt, dad die Konstanten <Ma- 
Irixidemenlel von unterschiedlieher GiblJenordnunR 
sind: das bewirki vine weitere ivesentlirhe Fehlererhd- 
hom;. ix-i der BererdtnunR der Konstanten. Auth dutch 
eine liildunRsvorschrifl fiir p, die von der Wahl der 


Zeitkonstanten der ObertraRunRsfunktion abhanRiR isl, 
isl das Problem nicht enlRulliR losbar. 

Aus die.sem Grunde wurde im RechnerproRramm cine 
interne ZeitnormierunR vorRenommen. wolx-i die Nor- 
mierunR.szeit vt>n den f:inRalx>daten abhiinRiR 

isl. Fiir p wird eine von den EinRabedaten unabhiinRiRc 
BildunRsvor.sctirifl bcnulzt. .So wird es miiRlich. die Elc- 
mente der Koeffizienlenmntrix dcs linetiren GleichunRs- 
sy.stems in rechcntechni.sch RiinstiRcn Bereichen zu hal- 
ten. 

4. Bemerkungen zur Programmorganisatlon 

Das ProRramm wurde in dcr m.a.-,chincnoi-ienlierten Pro- 
Rrammiersprache fiir den Klcinrechner C 82(1:5 Re- 
.schrieben. WeRen des RiolJen Speicherplatzbedarfs 
wurde eine AufieilunR in zwei ProRrammicile vorRe- 
nommen. 

Zu BeRinn dcs 1. Programmteil.s muC dcr Rcchner 
mit folRcnden UnlerproRrammen Reladen werden: n-le 
Wurzel aus x. c'’-ProRramm, In x-Programm. linearcs 
Gleichung.s-sy.stem. Danach wird der Lo(iislreifen fur den 
Prngrammablauf des I. Programmteil.s {Hauptpro- 
gramm), dann der Einualjestreifen mit den aktuelien 
Parametern einRclc.sen. Nach der BeendiRunR dcs ersten 
ProRrammleils wird ein Eochstreifen mit den im Unter- 
proRramm fiir die LiisunR dcs linearen Gleichungs- 
systems erreehnclcn Konstanten in interncr DarstellunR 
ausRCReben. 

Fiir den zweilen ProRrammteil werden folgcnde Un- 
terproRramme benbtiRl: Wurzel aus x. e^-ProRramm, 
In x-Pi'oRramm, sin x. cos x, arct.in x..Nach der 
dung des Retainers mil die.sen UnlerproRrammen wird 
das HauplproRramm dcs 2. ProRrammleils und der am 
Endc des 1. ProRrammleils nusRedruckto Lochslreifen 
eingelesen, die AusRabe erfoiRt nach BeendiRunR des 

2. ProRrammleils iiber das .Sehreibwerk des Rechners. 
Als aktuelle Parameler werden einRcRcben: 

I. Pammefer der UberlraRunRsfiinklion im Ziihler und 
im Nenner. Vielfachheilcn der Pole. .SehwinRunus- 
Rlicder usw. 

Gewiinschle Systemreaktion (Olx-rRanRsfunktinn. 
Gewichtsfunktion o. :i.) 

3. Gewiinschtc DarslellunRszcit der Syslemronktion und 
Sdirillweite. 

.5. Ergebn.sse 

Das Verfah.en wurde an einer Viclzahl der verschie- 
clensten Beispielc Retestel. Dazu wurde zur OberprilfunR 
die Riiektransformation zu.siitzlich noch auf dem meisl 
iiblichen W'eRc mit Korrespondenztabelle durchRefiihrt. 
f>ie so erhallene Zeilfunklion wird vcrcinf.ichend .-ils 
die exaktc Zeitfunklior. bezeichnet. deren Wertcverlauf 
mitlels dcr vorhandenen UnIcrproRramme ebenfalls 
iiber den Rechner ermittclt werden kann. 

Die (lurch das ProRramm crrechneten Zeitfunktionen 
wiesen bei alien Beispielen nur .sehr RcrinRo Abwei- 
chunRen auf. Der ZusammenhanR zwischen waehsendem 
Fehler und waehsendem Nennergrad von X., (p) wurde 
besliitlRt. Dennoch trat bei einem Nennercrad von sicben 
ein Fehler. bezogen auf den Rrubten Funklionswerl, von 
nur 3 ■ lO"' auf. Bei einem Nennergrad von elf war das 
Pn'Rramm ullerdinRS uberfordert. die ausRedruckle 

Wertetabelle lieferte nicht einmal Anhallspunkte uber 
den reaie.T t uniuiutisvenaui vuii x^.i;. 




Aus dem Vergleich des vorgestellten Algorithmus mit 
anderen, aus der Lfiieratur ut^Kaiiiinii vciiamt-'ii i->— 
geht hervor. daO dieses Programm cinen fiir technische 
Anwendungen gunstigen KnmpromiO zwischen erfor- 
derlidu-n Cenauigkeil und rechentechnischem Aufwand 
darslellt. 


Zasammenfassung 

In dieser Arbeit wurde ein Programm^ur numerischen 
Laplace-ROii^ransformation vorgcstglU. Der angegebene 
Algorithmus wsruht auf dem Vcmhrcn der Zerlegung 
der Bildtunklioirsin Partialbndoie, die dann getrennt 
• iicktransformiert wstrden. Summe der entslehenden 
Teilfunktionen ist diesesuimte Zeitfunklion. 

Auftretende Problem^^ervorgerufen durch redien- 
technisch ungiinstige rarmebs^eines zu idsenden Clei- 
chungssystems. werdm durch ' entsprechende Normie- 
rung mit einem WIfsparameler weilgehend einge- 
schrankt. Die pralaische Anwendungsgrenze liegt bei 
einem Ncnnergrad/der BildfunkUon von sioben bis acht. 

PesMMe 

B pa6ore npcACTaBaciia iiporpaMMa iiiicbpoBOH oOpaTHoii 
TpanabopManmi no Ranaacy. Raii^iii aaropHTM ocHOsaH 
na pa3ao/KCHMn''>li4;nKumi OTO(^)KCHMa iia ajiCMCiiTap-, 
iiHc Apof)n. Kawflas Nn kotojwx nocac aroro OTgcabiio 
dOpariio Tpanc(bopMiipyb<^c«.y€yMMa o6pa30BanHi.ix nog- 
cjtyiiKiim'i saaxeTCx iicKO.MOj/f{)y'"‘'tMcri BpCMCiin. 

npoOacMM. BOamiKaKMuyc ih^-aa TpygiiocTM nenoAbaoBa- 
mix BbimicAiiTCALiiori TCxmiKii, B^bisacMOM iicdaaronpHXT- 
iibiMii ()}op.MaMit pcmnc/ioii oictcmu ypaBiiCHnii. aiiaMti- 
TCAbllO yMCIIblliaiOTCfl CniO.MOlUbK) COOTBCTCTBCHHOl'l liOpMH- 
PORKK CO BCnOMOraTcAllblM Iiapa.MCTpOM. npaKTHMCCKllfl 
iipcgCA iipiiMCHcmix HajmgiiTCB npn crcncmi airaMctiaTCAB 
d>yiiKiuni OTOf>pa>KCHHa or 7 go 8. 


AbschlieDend mwiiten die Verfasser auf dicsem Wege 


Sektion Schiffstcchnik fiir ihre Milarbcit bei der F’ro- 


grnmmierung und fiir die DurchfUhrung der numeri- 
schen Redinungen im Lchrreciincrkabinetl der Sektion 
Schiffstcchnik danken. 
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Summary 

This paper presents a program for the numerical 
Laplace return transformation. The algorithm described 
is based on the methdd of breaking the picture function 
down to form partial fractions whicli arc then trans- 
formed back separately.. The sum of the resulting par- 
tial functions is the sought time function. 

Any problems caused by the forms of a .system of 
equations being unfavourable due to the computer con- 
figuration are restricted to a large extent by appropriate 
standardisation using an auxiliary parameter. The prac- 
tical limits for use of the program are reached when 
the degree of the denominator of the picture function 
lies within the range of seven to eight. 

Resume 

L’article pr6senfe unXprogramme de^fransformation 
nun Unique inverse de DHplace. L’algorithme donn6 sc 
ba.se sur le proc6d6 de cWompogUion de la fonclinn 
d’irr.age en fractions particlliV ti^i .sont ensuite relrans- 
form^cs s6par6ment. V 

La somme des fonctions Mrtbdlos obtenucs est la 
fonction de temps cherch6e. / N. 

Des problcmes cau.sds par des Rirrntss do calcul d6s- 
avanta'geuses d’un syslernfc d'dqunlion uN:6.soudre .sont 
largement rfeduits gr5ce/a la standardisation adequate, 
effectude a I’aidc d'un/paramcfre auxiliaire. La limiie 
d’applieation pratique fe .silue a un degre de denomina- 
teur de la fonction d'image enire sept et huit. 
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